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spectral  density  of  energy  distribution  of  turbulence  motion 
or  of  the  energy  spectra  of  turbulence  as  formulated  by  Obukhov-  i 

Kolmogorov.  He  also  reviews  various  experimental  studies  con¬ 
firming  the  presence  of  buoyancy  forces  in  the  spectrum  of  the  { 

turbulence  interval  of  wave  numbers.  Experimental  curves  ar.d  j 

equations  are  presented  showing  that  buoyancy  subregions  exist 
in  the  spectrum  of  developed  continuous  turbulence  during  stable 

stratification  of  the  atmosphere.  In  the  case  of  stable  wave 
disturbances  the  spatial  energy  flux  retains  its  direction  towards  i 

large  wave  numbers.  The  energy  arriving  in  large  scales  in  trans-  • 
mitted  to  the  region  of  submolecular  scales.  In  this  case  a  ' 

discrete  spectrum  is  obtained.  If  the  medium  possesses  viscosity 
a  "pumping  out"  of  energy  occurs  from  any  mechanical  notion  • 

forming  in  this  region.  In  the  case  of  very  strong  thermal 
stability  of  the  medium,  the  spectrum  of  turbulence  may  he 
lzed  in  a  wide  zone  of  wave  numbers.  In  the  case  of  less  tr.-.rma. 
atmospheric  stability,  as  in  the  case  of  large  dynamic  Instability, 
and  also  in  the  region  of  wave  numbers  far  from  the  subregion  of 
buoyancy,  the  turbulence  spectrum  obtained  is  considerably  wider 
and  differs  substantially  from  spectra  of  purely  harmonic  fluctua¬ 
tions.  The  most  frequent  case  is  tout  of  a  continuous  spectrum 
of  turbulence  in  which  kinetic  energy  is  transferred  from  smaU  » 
to  large  scales.  Prig,  art,  has:  5  figures,  11  formulas. _ I 
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(U)  The  results  of  aircraft  investigations  of  atmospheric  turbulence  • 

between  IOC  ra  to  3  Km  in  the  area  around  the  city  of  Dolgoprudnyy 
carried  out  by  the  Central  Aerological  Observatory  in  1946  are 
analyzed,  A  meteorograph  and  an  accelerograph  for  recording  accele¬ 
ration  in  a  vertical  plane  were  installed  on  the  aircraft.  During 
the  flights  vertical  soundings  were  made  from  areas  at  heights  of 
100,  300,  500,  1000,  1500,  2000,  and  ?G0Q  m  and  at  the  surface  of 
the  earth.  The  turbulent  exchange  (K)  was  computed.  Individual 
flights  analyzed  with  the  aid  of  graphs  of  the  vertical  profile  of 
the  turbulence  coefficient  and  the  daily  variation  of  the  turbulence 
coefficient  at  different  altitudes.  Under  anticyclonic  conditions 
the  turbulence  coefficient  has  a  distinct  vertical  profile  and  its 
maximum  is  around  13  hours  at  all  levels.  Frequently  one  to  two 
maxima  were  observed  and  their  values,  are  determined  by  the  thermal 
stratification  of  the  atmosphere.  Trie  turbulence  coefficient 
manifests  a  daily  variation.  The  maximum  K  values  are  reached 
simultaneously  at  all  levels  in  the  cay time,  but  the  beginning 
of  daily  variations  is  most  clearly  manifested  at  100  m. 
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ABSTRACT 

(U)  The  deductive  turbulence  theory  proposed  by  Chandrasekhar  is 
based  on  postulates  of  homogeneity,  isotropy  and  stationary  state 
and  also  on  the  hypothesis  of  quasi-normality  (hypothesis  of 
Millionshchikov)  associating  correlation  tensors  of  the  second 
and  fourth  rank.  In  this  paper  the  authors  generalize  the  equation 
of  Chandrasekhar  for  the  case  of  non- stationary  turbulence.  There 
is  derived  a  dynamic  equation  for  determining  the  scalar  correla¬ 
tion  tensors  describing  a  two  point  space-time  c  'relation  of  the 
velocities  of  turbulent  flow  and  constituting  a  generalization  of 
the  equation  of  Chandrasekhar  for  stationary  turbulence.  In  the 
derivation  there  were  assumed  conditions  of  homogeneity  and  isotrop 
and  also  the  generalized  hypothesis  of  cuasi-ncrmal  distribution  of 
velocities.  In  the  case  of  the  direct  derivation  of  the  stationary 
equation  there  is  obtained  a  pair  of  equations  that  have  not  teen 
considered  by  Chandrasekhar,  But  these  equations  have  to  be  disre¬ 
garded  because  they  cannot  be  obtained  from  a  non-s tat ionary  equa¬ 
tion.  In  principle  it  is  pcssibie  to  obtain  another  equation  for 
non- stationary  turbulence  which  has  not  yet  been  investigated.  The. 
generalized  equation  of  Chandrasekhar  correctly  describes  the  last 
stage  of  the  origin  of  turbulence.  The  stationary  equation  or 
Chandrasekhar  for  the  soatial- temporal  longitudinal  correlation 
function  can  be  solved  by  representing  the 'solution  in  a  power 
series,  Orig.  art.  has:  2  figures,  112  formulas, 
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DAILY  VARIATION  IN  THE  TURBULENCE 
COEFFICIENT  ABOVE  FLAT  TERRAIN 


V.  D.  Litvinova 


Results  are  given  of  aircraft  investigations 
of  turbulence  around  Dolgoprudnyy  from  an  altitude 
of  100  m  to  3  km.  It  was  found  that  under  anti- 
cyclonic  conditions  the  maximum  value  of  the 
turbulence  coefficient  at  all  levels,  is  observed 
around  1300  hours .  The  value  of  the  maximum 
depends  on  the  thermal  stratification  of 
atmosphere.  Diurnal  variations  in  the  turbulence 
coefficient  are  most  noticeable  at  a  height  of 
100  m. 


For  investigation  of  a  number  of  physical  processes  in  the 
atmosphere  it  is  necessary  to  know  the  intensity  of  turbulent 
exchange  of  heat,  moisture,  and  momentum. 


Experimental  investigations  of  the  conditions  of  appearance 
and  distribution  of  turbulence  were  conducted  at  TsAO  (Central 
Aerological  Observatory),  GGO  (Main  Geophysical  Observatory), 
and  others.  For  free  atmosphere  the  most  detailed  investigations 
of  variability  of  the  turbulence  coefficient  with  time,  height, 
and  horizontal  extent  are  given  in  the  works  of  P.  A.  Vorontsov  fl] 
N .  7.  Pir.us  [5],  and  V.  Ye,  Minervin  [4], 


In  this  work,  according  to  data  from  flights  conducted  by 
TsAO  in  1946  on  a  PO-2  aircraft  in  the  Dolgoprudnyy  area,  certain 
characteristics  of  turbulent  exchange  in  the  lower  troposphere  are 
examined . 
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Preceding  Page  Blank 


The  aircraft  carried  a  meteorograph  ana  an  accelerograph  of 
type  SP-11,  intended  for  recording  g  of  th?  aircraft  in  the  vertical 
plane.  During  flights  vertical  sounding  was  made  of  the  atmosphere 
with  platforms  at  heights  of  100,  300,  500,  1000,  1500.  2000,  and 
5000  Tt  above  the  earth's  surface.  The  durar.ion  of  recording  of 
accelerograpb  on  these  platforms  did  nor.  exceed  one  to  two  minutes. 
In  ail  48  flights  were  made,  one  or  two  flights  a  day,  with  four 
series  of  quickened  soundings  of  the  atmosphere  being  made  on  three 
to  five  flights  daily. 


In  1965  all  data  from  these  flights  was  reanalyzed,  coefficients 
of  turbulent  exchange  were  calculated,  and  graphs  of  dependence 
of  the  coefficient  of  exchange  on  heights  and  time  of  day  were 
plotted.  The  coefficient  of  turbulent  exchange  was  calculated  by 
the  Lyapin-Dubov  formula,  definitized  tc  M,  A.  German  in  work  [2*: 

„  VJWiAt* 

A - - 

where  V  —  average  flight  speed  of  aircraft  on  given  platform, 
m  —  average  absolute  velocity  of  vertical  gust  of  air,  during 
the  calculation  of  which  the  mechanical  and  flight  characteristics 
of  the  aircraft  were  considered,  &  —  average  time  of  preservation 
of  sign  of  pulsations  in  speed,  t—  value  of  transfer  function  for 
the  given  type  of  aircraft,  given  in  work  [5]. 


Below  will  be.  examined  only  cases  of  quickened  sounding 
of  the  atmosphere. 


Fli ghts  of  11  August  194 6 


On  11  August  1945  three  f Sights  were  made  at  13~7,  1552,  ana 
1830  hours.  Flights  passed  through  the  centra]  part  of  an  anticyclone 
Variable  overcast  with  weak  north  and  nor:.:, west  winds  was  observed. 
Platforms  were  at  heights  of  100,  500,  1000,  .00<\  and  >000  m. 


Pig.  1.  Vertical  profiles  of  the 
turbulence  coefficient,  a)  II  Aug  1946, 
b'l  12  Aug  1946,  c)  29  Sep  1946, 
a)  25  Oct  1946. 

As  can  be  seen  from  Fig.  la  during  the  flight  of  the  aircraft 

^  2 

at  1317  hours  the  value  of  K  at  a  height  of  100  m  was  24  m  /r ,  at 
500  m  —  18  m2/s,  and  at  1000  m  —  21  m2/s.  The  value  of  K  then 
decreases  rapidly  to  2  rr//s  at  a  height  of  2  km.  At  a  height  of 
3  km  the  flight  was  calm.  The  vertical  temperature  gradient  varied 
from  0.7  to  0.9°/100  m.  At  a  height  of  2.8  km  there  was  an  inversion. 
A+  &  height  of  1,5  km  was  the  lower  FrCu  and  Cu  cong  boundary. 
Unfortunately,  no  g  recording  was  made  at  this  height. 

In  155'2  hours  the  value  of  the  turbulence  coefficient  at 
2  2 

100  m  was  18  m  /s,  at  500  m  —  15  m  /c,  and  from  1000  m  the  flight 
was  calm.  The  vertical  temperature  gradient  to  a  height  of  66 0  m 
was  1°/100  m,  and  at  higher  altitudes  0.75°/100  m.  At  2.8  km 
7  =  0,1/100  m. 

At  1830  hours  the  turbulence  coefficient  at  100  n.  was  12  m  /s, 
and  at  500  m  the  flight  was  calm.  The  value  of  7  changed  from 
0.75  to  1°/100  m,  and  at  2.8  km  a  temperature  inversion  was  observed . 

As  can  be  seen  from  Fig.  la,  maximum  intensity  of  turbulence 
was  observed  at  1300  hours,  and  it  gradually  weakened  toward 
l8C0;  at  a  height  of  ICO  m  the  value  of  the  turbulence  coefficient 
decreased  frcm  24  m"/s  (at  1300'  to  12  m2/ s  (at  1800),  anu  at 
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^00  m  decay  of  turbulence  toward  evening  v/sg  still  faster:  from 
23  m^/s  at  1300  to  15  m2/s  at  1500,  and  at  1800  at  this  height 
flight  was  absolutely  calm.  As  can  be  seen,  changes  of  turbulence 
coefficient  with  height  follow  well  the  changes  with  height  of  the 
vertiial  temperature  gradient,  especially  in  daytime  houisc  With 
decrease  of  the  vertical  temperature  gradient  of  the  value  cf  K 
decreases . 


2 •  Flights  of  12  August  1946 

Flights  passed  through  the  central  part  of  an  anticyclone 

during  scantily  clouded  weather  with  weak  northwest  winds  at  the 

surface.  Flights  were  made  at  0552,  0848,  1139,  and  1736  hours  at 

heights  of  100,  5^0,  1000,  2000,  and  3000  m.  The  turbulence 

coefficient  was  calculated  for  only  the  last  three  flights,  since 

the  first  was  calm.  At  100  m  at  0848  hours  the  value  of 
2 

K  =  .16  m  /s  (Fig.  lb).  In  the  layer  to  300  m  were  observed  raised 
values  cf  vertical  temperature  and  wind  gradients  (7  =  0.92°/lC0  m, 

2  =  0r83  m/s/100  m) .  In  the  300-2560  m  layor  values  cf  7  did  not 
exceed  0.85°/100  m,  and  3  was  not  over  0.35  m/s/100  m.  At  2560  m 
was  an  inversion. 

Further  flight  in  the  entire  layer  sounded  was  calm.  In 
1139  hours  the  on  vertical  profile  of  the  turbulence  coefficient 
two  maxima  are  observed.  The  first  maximum,  caused  by  the 
influence  uf  the  underlying  surface,  is  at  a  height  of  100  mj  the 
value  cf  K  here  is  16  m  /s .  The  value  of  the  vertical  temperature 
gradient  is  0.86 '""/ICO  a,  and  that  of  wind  is  7.1  m/s/ioo  n.  At 

500  m  the  value  of  K  is  11  m  /s,  and  then  again  increases,  attaining 

2  0 
18  r.  /s  at  a  height  of  100  m.  Appearance  of  the  second  maximum  of 

the  turbulence  coefficient  can  be  explained,  apparently,  by  thermal 
stratification  of  the  atmosphere.  For  this  layer  there  characteris¬ 
tically  is  presence  of  a  large  value  of  vortical  temperature  gradient 
(7  =  C.97°/100  m)  and  sharp  decrease  of  vertical  wind  gradient 
(to  0.35  m/s/100  in).  Above  this  level  flight  was  calm.  By  1736  hours 
there  is  decay  of  turbulence  in  the  entire  sounded  layer.  Values 
of  K  at  all  heights  to  1000  m  by  this  time  no  longer  exceeded 
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several  m/s.  At  heights  of  2000  and  3000  si  flight  was  calm.  The 
vertical  temperature  gradient  in  the  300-1300  m  layer  equaled 
0.82-0. 94 °/100  m  and  at  the  surface  of  earth  and  at  high  altitudes 
was  0.64-0.73°/100  m.  The  value  of  £  to  a  height  of  3000  m  was 
2.7 2  m/s/100  m,  and  further  0.1-0. 3  m/s/100  m. 

3 .  Flights  of  29  September  1946 

Flights  passed  through  the  spur  of  an  anticyclone.  Cloudless 
weather  with  weak  north  winds  at  the  surface  was  observed. 

Data  obtained  in  this  day  are  the  most  interesting.  During 
the  day  six  flights  were  made.  Therefore  it  was  possible  to 
trace  most  fully  the  atmospheric  composition  to  a  height  of  2  km. 
Flighs  were  made  at  0640,  0338,  10^9,  1238,  1433,  and  l62t  hours 
at  heights  of  100,  300,  300,  1000,  1500,  and  2000  m. 

According  to  data  from  five  flights  vertical  profiles  of  the 

turbulence  coefficient  were  plotted.  For  the  1625  flight  the 

turbulence  coefficient  was  not  calculated,  since  turbulence  was 

*  «• 

not  observed  continuously,  but  only  at  times. 

A:s  can  be  seen  from  Fig.  1c,  at  0640  the  turbulence  coefficient 

at  100  m  was  17  m  /s,  while  at  higher  altitudes  flight  was  calm. 

In  the  layer  of  atmosphere  from  the  earth  to  400  m  a  deep  inversion 

was  observed.  By  0838  the  temperature  inversion  had  vanished  and 

the  vertical  temperature  gradient  increased  rapidly  to  1.00°/100  m, 

while  the  wind  gradient  increased  to  3.75  ir/s/lOO  m.  At  100  m 

X  =  12  irZ/s.  By  1029  the  vertical  profile  of  K  significantly 

changes.  On  the  figure  is  seen  a  well-defined  maximum  of  the 

turbulence  coefficient  in  the  layer  from  the  earth  to  300  m,  where 
2 

K  =  14-19  m  /s.  At  heights  from  1000  m  to  £000  m  flight  was  calm. 
Values  of  vertical  temperature  gradient  to  a  height  of  1100  m 
varied  from  0.8  to  1.0°/100  n,  and  wind  gradient,  varied  from 
0.2  to  1.0  m/s/ 100  m.  By  1238  hours  atmospheric  turbulence  was 
strengthened  sharply.  On  the  profile  of  the  turbulence  coefficient 


2 

two  maxima  are  observed.  At  a  height  of  100  m  K  =  4o  m  /s,  and 

then  there  is  a  certain  tendency  toward  slight  decrease  at  a 

height  of  300  m,  after  which  the  value  of  K  increases  to  500  m, 

2  2 

where  it  attains  35  m  /s .  At  1000  m  K  ■  31  m  /s .  The  vertical 
temperature  gradient  in  this  layer  is  equal  to  1°/1000  m,  and 
that  of  wind  is  0.15-0.58  m/s/100  m.  At  heights  of  1500  and  2000  m 
flight  was  calm. 

By  1433  the  vertical  profile  of  K  has  more  smoothed  form.  In 
the  layer  from  100  to  500  m  the  turbulence  coefficient  was 
20-24  nF/fe  and  at  1000  m  was  10  m2/s .  In  this  whole  layer  the 
vertical  temperature  gradient  is  equal  to  1°/100  m,  that  of  wind 
to  0.15-1.00  m/s/100  m.  At  2000  m  flight  was  calm.  By  1625  almost 
full  decay  of  turbulence  had  taken  place. 


Fig.  2.  Daily  variation 
of  turbulence  coefficient 
for  various  heights. 


In  Fig.  2  is  represented  daily  variation  of  turbulence 
coefficient  at  heights  of  100,  500,  and  1000  m.  As  can  be  seer, 
from  the  figure,  by  0900  hours  turbulence  extends  to  an  1  r.n ! gnl  f  leant 
height,  with  values  of  K  »  12-17  m2/s  being  attained  at  a  nelght 
of  100  m.  When  at  100  m  turbulence  is  developed.  It  starts  to* 
spread  upwards.  Although  turbulence  at  100  rr.  starts  to  develop 
considerably  earlier  (around  0600),  it  attains  maximum  at  all 
heights  simultaneously  from  1200  to  1300  hours. 
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The  value  of  the  turbulence  coefficient  during  this  rime  is 

O 

li o  m^/s  at  IOC  m,  gradually  dropping  with  height  to  J1  m  /s  at  1000  m. 


Then  the  intensity  of  turbulence  decreases  and  at  1433  K  changes 
from  24  m^/s  at  H  =  100  m  to  10  m^/s  at  H  =  1000  m. 


For  more  detailed  consideration  of  this  effect  in  Fig.  3 
it  is  shown  how  the  value  of  K  changes  with  time  at  various  heights 
in  daylight  hours. 


Fig,  3-  Time  cross  section  of  the 
turbulence  coefficient  for  29  Sep  1946. 


Unfortunately,  there  are  no  night  observations.  Apparently, 
at  nignt,  when  in  the  lower  part  of  the  troposphere  temperature 
stratification  becomes  stabler,  the  intensity  of  turbulence  should 
decrease  sharply. 

4 .  Flights  of  25  October  1946 

Flights  passed  through  the  rear  part  of  a  cyclone,  where 
secondary  fronts,  displaced  to  the  south  were  observed.  In 
the  Moscow  region  were  observed  stratus  clouds,  weak  snowfall  at 
times,  and  northwest  wind . 


this  day  four  flights  were  Trade :  at  092'., 


1336  hours.  Platforms  were  made  at  lower  and  upper  cloud  boundaries. 


i 

i . 


At  0740  at  heights  of  iOO  and  I/O  m  (at  trie  -ower  cloud  boundary) 

the  aerologist:  noted  moderate  bumping  of  the  PO-2  aircraft.  The 

2  2 

turbulence  coefficient  was  12  m  /s  at  100  m  and  15  m  /o  at  250  m 
(Fig.  Id).  At  remaining  heights  flight  was  calm. 

By  1128  turbulence  had  strengthened.  The  turbulence  coefficient 
was  considerable  from  the  lower  flight  platform  to  the  height  of 

p 

the  lower  St  boundary,  which  was  at  450  m.  At  100  m  K  was  16  m  /s, 

O 

while  at  300  m  the  value  decreased  to  11  m  /s,  and  then  increased 

p 

at  the  lower  St  boundary  to  18  m  /s .  The  vertical  temperature 
gradient  in  this  layer  varied  from  0,33-0.78/100  m.  The  aerologist 
also  noted  moderate  bumping  at  three  lower  levels.  Strengthening 
of  intensity  of  turbulence  at  the  lower  cloud  boundary  agrees  well 
with  data  from  work  [4]. 

At  1335  overloads  of  the  aircraft  were  observed  only  in  the 
form  of  bursts,  and  it  was  not  possible  tc  calculate  the  value 
of  the  turbulence  coefficient. 

As  a  result  of  the  foregoing  the  following  preliminary 
conclusions  can  be  made. 

The  turbulence  coefficient  under  anticycionic  conditions  has 
well-defined  vertical  profile.  On  it  most  frequently  are  noticeable 
one  or  two  maxima,  the  values  of  which  are  determined  by  conditions 
of  thermal  stratification  of  the  atmosphere . 

Values  of  the  turbulence  coefficient  have  daily  variation. 
Maximum  of  K  is  attained  in  the  daytime  hours  at  all  heights 
simultaneously,  but  the  beginning  of  diurnal  variations  is  most 
clearly  expressed  for  the  height  of  100  m. 

Bibliography 


1.  Vorontsov  ?.  A.  Aerologicheskiye  issiedovar.iya  pogranichnogo 
sloya  atmosfery.  Gidrometeoizdat  (Aerciogicai  investigations  of  the 
atmospheric  boundary  layer).  L.,  19 60. 


■o 


8 


2.  German  M.  A.  0  turbulentnom  obmene  v  oblakakh  (On  turbulent 
exchange  in  clouds).  Meteorologiya  i  Gidrologiya,  No.  10.  M.,  1965. 

5.  Dubov  A,  S.  Opredeleniye  koeff itsiyenta  turbulentnogo 
obmena  po  uskoreniyu  sajjoleta  (Determination  of  the  coefficient  of 
turbulent  exchange  in  accordance  with  acceleration  of  an  aircraft). 
Trudy  GGO,  vyp.  98,  M,  1959. 

A.  Minervin  V.  Ye.  Turbulentnost 1  v  nizhnem  sloye  atmosfery 
i  oblakakh  nizhnego  yarusa  (Turbulence  in  the  lower  layer  of 
atmosphere  and  clouds  of  the  lower  stratum).  Trudy  TsAO,  vyp.  Jl. 

M . ,  1966 . 

5.  Pinus  N.  2.  Nekotoryye  osobennosti  razvitiya  turbulentnosti 
naa  ravninnoy  mestnost'yu  (Certain  peculiarities  of  development 
of  turbulence  above  flat  terrain).  Izvestiya  AN  SSSR.  Fizika 
atmosfery  i  okear.a,  t.  1,  No.  3,  M.,  1965. 


Submitted 
15  August  1966 


\ 


\  ! 
\  i 
i  t 

j 

.  1 


9 


THEORY  OP  SPATIAL-TEMrORAL  CORRELATION 
OF  VELOCITIES  IN  AN  ISOTROPIC 
TURBULENT  FLOW 


V.  I.  Smirnov  and  B.  Sh.  Shapiro 


Conditions  of  homogeneity  and  isotropy,  and 
also  the  hypothesis  on  quasi-normality  of  velocity 
distribution  are  used  for  derivation  of  an  equation 
for  the  Eulerian  two-point  space-time  correlation 
function  in  nonstationary  turbulent  flow,  general¬ 
izing  equation  the  Chandrasekhar  for  stationary 
turbulence.  Solution  of  the  derived  equation  for 
the  last  stage  of  degeneration  of  turbulence  agrees 
with  results  obtained  by  other  authors.  A  method 
of  solving  the  Chandrasekhar  equation  by  expanding 
the  solution  in  power  series  is  developed. 
Coefficients  of  the  series  are  determined  uniquely 
if  space  and  time  correlation  functions  are  given. 
With  rapid  decrease  of  derivatives  of  the  corre¬ 
lation  function  with  distance,  it  is  sufficient 
to  assign  the  space  correlation  function.  The 
time  microscale  of  correlation  is  determined  by 
the  mean  square  of  velocity,  kinematic  viscosity, 
and  the  space  correlation  function.  There  is  a 
class  of  universal  solutions,  i.e.,  solutions 
not  depending  on  viscosity  and  intensity  of 
turbulence.  Certain  solutions  not  depending  on 
the  time  interval  are  investigated. 


1.  Introduction 

Turbulent  flow  might  be  simultaneously  uniform,  isotropic, 
and  stationary  (in  a  statistical  sense)  only  in  the  presence  within 
it  of  corresponding  sources  of  energy,  replenishing  the  decrease 
of  energy  of  turbulence  due  to  its  dissipation.  The  theory  of 
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turbulence  is  most  simply  built  around  the  three  above-mentioned 
conditions,  and,  bypassing  the  question  of  sources  of  energy, 
around  introduction  of  certain  more  or  less  formal  assumptions. 

In  1955-1956  Chandrasekhar  offered  the  deductive  theory  of 
turbulence  [9>  10],  based  on  postulates  of  homogeneity,  isotropy, 
and  stationarity,  and  also  on  the  well-knovm  hypothesis  of 
quasi-normality  (Millionshchikov  hypothesis),  associating  correlation 
'tensors  of  second  and  fourth  ranks.  On  the  basis  of  the  Navier- 
Stokes  -equation  for  incompressible  liquid  Chandrasekhar  derived  a 
nonlinear  partial  differential  equation,  containing  only  one  unknown 
scalar  function,  depending  on  space  and  time  intervals.  In  view 
of  the  assumed  difficulty  of  solving  this  equation  in  the  general 
case.  Its  author  limited  himself  to  finding  a  solution  for  the  case 
when  viscosity  does  not  play  a  role  and  can  be  considered  equal 
to  zero,  and  arrived  at  the  conclusion  that  theory  does  not  contradict 
the  statements  of  the  well-knovm  Kolmogorov  theory  on  the  inertial 
interval. 

In  developing  the  Chandrasekhar  method,  S.  Panchev  derived 
a  dynamic  equation  for  the  space-time  correlation  function  and  an 
equation  for  other  correlation  functions  [11], 

Later  A.  I.  Ivanovskiy  and  I.  P.  Mazin  showed  [2]  that  with 
application  of  the  above-mentioned  assumptions  one  more  equation, 
different  from  the  Chandrasekhar  equation  can  be  obtained.  The 
form  of  the  equation  derived  by  Ivanov  and  Mazin  was  definitized 
in  the  work  of  V.  I.  Smirnov  [5]. 

A  check  of  the  derivation  of  the  Chandrasekhar  equation, 
expounded  in  work  [9],  shows  that  in  compatations  error  in  sign 
is  a; lowed,  after  correction  of  which  instead  of  the  Chandrasekhar 
equation  a  pair  of  new  equations  is  obtained.  However,  by  modifying 
’the  derivation  one  can  also  obtain  the  Cr.andrasekhar  equation. 

It  is  shown  below  that  rejection  of  the  condition  of  stationarity 
makes  it  possible  to  prove  that  this  pair  of  new  equations  has 
no  meaning. 
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In  this  work  the  Chandrasekhar  equation  is  generalised  for  the 
case  of  nonstationary  turbulence,  solution  of  this  equation  for 
the  last  stage  of  degeneration  of  turbulence  is  shown,  and  certain 
solutions  of  the  Chandrasekhar  equation  (for  the  case  of  stationary 
turbulence)  are  obtained. 

2 .  Equation  Describing  Nonstationary  Turbulence 

Rejecting  the  condition  of  stationarity,  adopted  by  Chandrasekhar 
[91 >  we  will  preserve  conditions  of  homogeneity,  isotropy,  and 
incompressibility  of  liquid.  Derivation  of  the  equation  will  be 
a  generalization  of  the  Chandrasekhar  derivation  [9].  Let  us 
introduce  the  following  correlation  tensors: 


f)Uj  (P,  Q  u„  (r* ,  f)  «•  uiu'jul , 

_ _  _ 

n-utujuiut .  (i) 

Here  r*  and  i*  —  radius  vectors  of  two  points  in  space;  t'  and  t"  — 
two  moments  of  time;  u^  —  component  of  velocity  vector  along  the 
i-th  axis  of  the  Cartesian  coordinate  system  ,  in  which  ui  =  0; 
p  —  pressure;  p  —  density.  Averaging  is  done  by  sets. 

Because  of  isotropy  and  homogeneity,  all  correlation  tensors 
depend  only  on  I  —  t"  —  r\  and  also  on  t'  and  t"  (on  —  *'|  in 

the  case  of  stationarity).  In  particular  Q,f  and  TIJJt  have  the 
following  form  [9]: 


+  2Q>*«. 

7'«..-7TVA-(i-r+3n(V/.+y»)+sl^v  (2) 

—  determining  scalars  of  tensors  Qy  ar.d  Tv  & 


where  Q  and  T 


depending  on  r,  t'  and  t:,$  8y  -  Kronecker  de3.ta;  primes  near  Q  and 
T  signify  differentiation  with  respect  to  r. 

The  basic  statistical  hypothesis  will  be  formulated  as 

.  Q//.#~Q»Qyi+  Q«Qa+  0y(0.  o,  f,  OQu(0,  0,  t\  Or  (3) 

which  is  a  further  generalization  of  the  hypothesis  adopted  by 
Chandrasekhar  in  work  [9]  and  generalizing  the  Millionshchikov 
hypothesis „ 


Let  us  consider  the  Navier-Stokes  equation  for  the  first  point 


(r\  n 


da',  d  d»' 


(“) 


We  multiply  (4)  by  and  average  for  the  set 


*)V 

d 

<*' 

~  Xm 

d 

d 

*/ 

X,  * 

Here  by  v*  is  understood 
scalars,  we  obtain  the  first  equation  for  scalars 


(5) 


(-£_,Z>.)Q_(ri  +  S)r.  <6> 


where  Dn  —  Laplacian  in  n-dimensional  space  in  case  of  spherical 
symmetry,  equal  to 


D, . 


d*  ,  K  ~  1  d 
dfi  +  r  ir' 


(7) 


In  order  to  find  the  second  equation  for  scalars,  vfo  will 
start  from  the  Navier-Stokes  equation  for  the  second  point  (?\  /*): 


*£ 

dT 


Ujut 
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We  multiply  (8)  by  u'tu’k  and  average.  We  then  find 


Xlkj’^t^iTatJ,  (9) 

where  we  designated 

(jo) 

Turning  in  (9)  to  determining  scalars,  we  obtain  the  second 
equation  for  scalars 


(11) 


where  X  ->  determining  scalar  of  tensor  XlktJ. 

As  Chandrasekhar  showed  [9],  from  the  hypothesis  of  form  (p) 
and  condition  of  isotropy  it  follows  that 


£(s+'i)*--2Q£c‘Q'  (12) 

Generalization  (5)  in  the  case  of  nonstationarity  does  not  hinder 
the  obtaining  of  equation  (12),  where  the  form  of  the  last  term  of 
the  right  side  of  (.3),  not  depending  on  coordinates,  does  not  play 
6  rcle,  inasmuch  as  in  derivation  of  (12)  we  encounter  only 
derivative  Q*,jt  with  respect  to  S,  see  (10). 

It  is  easy  to  check  for  the  existence  of  identity 


Applying  operator  to  (6)  and  considering  (15),  we 

find 


Hi 


(14) 


(v-^Xi->°')'J“-(ri+5>Y’ 

whence  with  the  help  of  (12)  we  finally  obtain 

We  must  emphasize  that  in  this  equation  it  is  impossible  to  consider 
t'  =  t",  which  follows  from  the  derivation,  see  transition  from  (4) 
to  (5)  ana  from  (8)  to  (9).  Owing  to  this  it  is  impossible  to 
make  the  transition  from  (15)  to  the  equation  for  Q(r,  t,  t), 
derived  by  Pancher  from  the  Xarrnan-Howard  equation  with  the  help 
of  the  quasi-normality  hypothesis  [11]. 

In  the  particular  case  of  stationary  turbulence  scalars  depend 
only  on  r  and  —  so  that  ^  Let  us  assume  for 

definitiveness  that  <*>/':  then  —  —  and  (15)  becomes  the 

ar  it  v  ' 

Chandrasekhar  equation  for  Q(r,  t): 

£-(-£-,,°9<J— (is) 


Thus  for  Q  we  obtain  a  unique  equation  (not  counting  the  equation 
obtained  by  the  other  method,  see  below),  when  under  the  assumption 
of  static  ?  rity  there  car.  bo  obtained,  besides  equation  (io),  an 
a<ld  il.ii  pair  of  equations 

(17) 


if  during  derivation  of  the  second  equation  for  scalars  we  start 
with  the  Navier-Stokes  equation  for  the  first  point  (r1,  t'). 
Incidentally,  it  is  exactly  in  this  way  that  Chandrasekhar  proceeds, 
but  he  allows  error  in  sign  (equality  (25)  in  work  [9])  and  therefore 
arrives  ct  equation  (16).  With  nonstationarity  it  is  impose-'  r,le  to 
start  both  t/imes  from  the  first  point,  since  we  then  on-air.,  in 


lb 


v 


particular,  different  correlation  tensers  and  u‘,u‘kttj  with 

different  determining  scalars.  Thus  it  is  impossible  here  to  obtain 
a  nonstationary  equation  wnich  would  become  (17)  in  the  case  of 
stationarity.  Consequently,  the  pa.ir  of  equations  (17)  is  excluded. 

If,  in  contrast  to  (4),  we  start  from  the  second  point,  we  again 
obtain  equation  (15)  for  the  nonstationary  case  and  equation  (16) 
for  the  stationary. 

As  was  shown  by  Ivanovskiy  and  Mazin  [2],  and  also  by  Smirnov  [5], 
assumptions  of  homogeneity,  isotropy,  and  stationarity  permit  obtaining 
one  more  equation,  different  from  equation  (16).  Again  rejecting  the 
condition  of  stationarity,  we  leave  on  the  left  in  equations  (4)  and 
(8)  only  the  partial  derivative  with  respect  to  time,  multiply  left 
and  right  sides  and  average: 

*  (18) 


where  f  ^  t“ ,  as  in  (15),  and  <!>  -  is  a  function  whose  expression 
through  Q  has  not  yet  been  found .  However,  in  the  particular  case 
of  strongly  degenerated  turbulence,  when  in  equations  (4)  and  (5) 
it  is  possible  to  reject  nonlinear  members,  the  form  of  «l>  is  easy 
to  find 


*Q 

■■■  «e* 

et’dl" 


>s0iQ. 


(19) 


Ivanovskiy  and  Mazin  consider  [2]  that  nonlinear  members  can 
be  rejected  and  for  stationary  turbulence  also,  if  we  limit  ourselves 
to  tne  viscous  subdomain,  where  Re<l.  We  then  obtain  [5j 


(£  +  .’o!)q-o,  <?0> 

which  differs  from  (16)  if  we  reject  the  nonlinear  right  part. 
Apparently,  requirements  of  homogeneity  and  isotropy,  together  with 
the  hypothesis  of  qua&i-normality,  are  not  so  strong  as  to  make  it 
possible  to  derive  a  unique  equation. 
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It  is  possible  to  doubt  the  legality  of  rejection  of  nonlinear 
members  during  derivation  of  equation  (20),  inasmuch  as  stationary, 
undamped  turbulence  exists  as  a  result  of  nonlinear  energy  transfer 
from  bigger  vortexes  to  smaller.  Therefore  it  would  be  desirable  to 
obtain  full  nonlinear  equation  (18)  and  to  study  its  solutions.  Until 
this  is  done  there  is  no  possibility  of  selecting  some  one  of  the 
two  equations  (15)  and  (18)  and  stating  that  it,  indeed,  describes 
the  structure  of  turbulent  flow. 

It  is  easily  confirmed  that  equations  (15)  and  (19)  correctly 
describe  the  final  period  of  degeneration  of  turbulence.  If  it  is 
considered  that  at  initial  moment  of  time 

Q  (r,  0,  0)  —  const  exp  j  ,  (21) 

Then  function 


satisfies,  which  is  checked  directly,  both  equation  (15)  with 
zero  in  the  right  part  and  equation  (19).  Solution  of  equation- (22) 
for  t  =  t’  »  t"  becomes  solution  Q(r,  t,  t),  found  by 
L.  G.  Loytsyanskiy  [5]  for  t— 00,  when  the  form  of  initial  condition 
(21)  does  not  play  a  role 


Q(r, /,  0  — const exp ,  (25) 

and  also  solution  of  the  Karraan- Howard  with  initial  condition  (21), 
found  by  Batchelor  and  Townsend  [8] . 


Below  will  be  examined  certain  solutions  of  the  Chandrasekhar  (l€). 
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y ,  Solution  of  the  Chandrasekhar  Equation  Depending 

on  Time  Interval 


Let  us  examine  equation  (16).  It  Is  convenient  to  introduce 
dimensionless  longitudinal  correlation  function  f(r,  t).  By 
definition 


(24) 


Substituting  (24)  in  (16),  we  obtain  an  equation  for  f 


►  («r~  '**)/-«* 


d 

dr 


(25) 


Function  f(r,  t)  should  be  an  even  function  of  r  and  t  owing 
to  conditions  of  isotropy  and  stationarity  [7).  Let  us  assume  that 
f(r,  t)  can  be  represented  in  the  form  of  a  power  series 


/(r.O-VTaS't*:  (26) 

Putting  this  series  in  equation  (16),  we  find 

£  ^ik(2k  -  2**  f  £;(*_  l)(*-2)(4i*  +  8<  +  3)fl,*X 

4-1  *-i  MM 

-  Xr,;  ;«V  2  V  +  (27) 

MM  i-S  *- 0 

Yfe  divide  this  equation  by  r,  make  linear  transformation  of 
indices  of  summation  in  such  a  way  that  summation  starts  each  time 
from  the  zero  index,  and  equate  coefficients  of  members  with  identical 
powers  of  r  and  t.  Then  we  obtain. 

a  - - S - |V(i  +  3)  (i  +  2)  X 

<  » 

X  +  l) <4i’  +  32 i  +  83 )a,+**  +  J  <P  +  2)  X 

^—0  0 

X  (p  +  D  (2p  +  7  *-f] ,  ( ?*J ) 
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Let  us  write  out  the  matrix  of  coefficients  a/t : 


at  <***•♦•  •  • 

®I0  fl»l  <*»»♦••  ♦  • 

<*30  ast  **»••«  Cj*  ♦  • 

•  ♦••••#•  • 

ait  •  • 


(29) 


According  to  (2:8),  coefficients  with  value  of  second  (temporal) 
index  >1  depend  on  a  finite  number  of  coefficients  with  lower  value 
of  second  index,  i.e.,  each  column  of  matrix  (29),  in  addition  to 
elements  of  the  first  line,  is  a  function  of  elements  of  coi.umns  to 
the  left  of  the  examined  column,  and  also  of  elements  of  the  first 
line.  Consequently,  if  we  assign  all  elements  of  form  a„,  a*,  (i,  k>0), 
then  by  formula  (28)  all  coefficients  are  uniquely  determined. 

In  other  words,  the  Chandrasekhar  equation  (25)  uniquely  determines 
the  space-time  correlation  function  f(r,  t),  if  space  f(r,  0) 
and  time  f(0,  t)  correlation  function  are  given. 

If  functions  f(r,  0),  f(0,  t)  are  assigned  graphically,  they 
can  be  approximated  by  power  series,  and  we  can  thereby  determine 
some  quantity  of  first  coefficients  If  K's  of  first 

coefficients  are  given;  2K  coefficients  a^,  are  necessary  in 
order  to  determine  all  ajk  to  the  left  of  straight  line 
in  (29). 

In  the  particular  case  of  k  =  0  (second  column  of  (29))  we 
find  from  (28): 


a.  . ,  -  (i  +  A)  (<  +  2){i  +  !) <4i’  +  32c  +  63)  X 

•  *  (<  -4-  1>L 

'  X  +  **  (  P  +  2)  (  p  4- 1)  (2p  +  7)  apk J>#  a,^  ,j , 


i.e,,  all  elements  of  form  an ,  with  the  exception  of  are  determined 
only  by  elements  of  form  om.  Where  i  =  0!we  have  the  simple  formula 


f  ’f 

[  \rt 


«tt  —  +  54 


(31) 


experimental  check  of  which  would  be  interesting. 

With  sufficiently  rapid  decrease  of  derivatives  f(r,  t) 
with  distance  for  unique  determination  of  f(r,  t)  one  longitudinal 
space  correlation  function  f(r,  0)  is  sufficient. 
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If  the  following  condition  is  met 


(■£■  -  fl| - 0. 


(52) 


which  seems  natural,  one  can  determine  all  coefficients  of  expansion 
of  f(fk  i),  a*,  i.e.,  the  function  f(0,  t).  Integrating  (25)  with 
respect  to  r,  we  find 


(-y-AS)/-  -»  [c<)-|/<r',  O-i-D./K,  Orf-'J.  (JJ) 


where 


cm- CM  cf. 


If  we  substitute  in  (33)  f(r,  t)  in  form  (26),  then, 
besides  (28),  we  obtain 

*****  2(*+I)(2*  +  l)  *  k>0'  • 


(34) 


(35) 


These  formulas  together  with  (28)  make  it  possible  to  success¬ 
ively  determine  all  k>Q*  Actually  «•,  equals 

fc.  (36) 


C,-f/(r.O)io,/(r,  0)*.  (37) 

For  it  is  easy  to  obtain  the  formula 

Ck— 4  S  J  £  £  <*  +  0  0  +  2)  (2i  +  7)at+Xfaltk_,r*+*+t  Ur.  (jQ) 

As  was  shovm  above.,  from  the  first  column  of  the  matrix  of 
coefficients  (29)  we  can  determine  the  second.  Further,  with  the 
help  of  (28)  from  elements  of  the  second  column  we  determine  all 
elements  of  the  third  but  The  a*  elements  we  find  from  (35), 

inasmuch  as  element  has  been  found,  while  C-^  can  be  found  from 
(38),  since  in  the  right  pert  of  (38)  are  only  elements  of  form 
r.  s  0.1,  which  have  already  been  determined.  Analogously  we  find 
elements  of  the  fourth  column  from  elements  of  the  third,  etc. 

Let  us  consider  formula  (36).  If  f(r,  0)  and  «*,  (36)  determines 
the  Eulerian  temporal  microscale  of  turbulence  *0,  equal  to 

(59) 

Since  as,  obviously,  we  must  have  Ot^O,  from  (36)  we  obtain  the 
lower  limit  of  energy  of  stationary  turbulent  flows 

(4o) 

Gi 

-Ixperimental  check  of  formulas  (36)  and  (40)  is,  from  our  point  of 
view,  of*  great  interest. 

Let  us  consider  the  example  when  f(r,  0)  has  the  form 


fir .  0)-« 


T 


(«) 


We  easily  find  that 


fe_  -  -  -  -w.  .— 


I 


,/ 


''  S  «w— -*pr.  «*— C#— -6a,*. 

s'  / 

Then  fr6m  (36)  and  (39)  we  obtain 


f«-^7 


AM**) 


(42) 


As  appraisals  founded  on  experimental  data  for  turbulence  with 
large  values  of  Re  show,  usually  [6] 


;  iS 


wini  )  "  7  fi 

0  9  (*f  9-RW  1 


(“3) 


where  e  is  rate  of  dissipation  of  energy  of  flow.  Here  we  use  the 
well-known  expression  for  /£  [7] 


(44) 


Then  from  (42)  we  obtain 


T  ■ 


(*5) 


Let  us  note  that  in  examined  example  formula  (40)  is  satisfied 
according  to  (43). 

It  is.  known  that  correlation  function  f(r,  0)  is  poorly  described 
by  a  formula  of  form  (41),  so  that  relationships  (42)  and  (45)  have 
to  be  very  approximate.  On  the  other  hand,  from  the  expression  for 
C0  (37)  it  follows  that  CQ  is  determined  basically  by  small  distances, 
where  f(r,  0)  and  its  derivatives  are  not  very  small.  Therefore, 
ther;-  Is  reason  to  think  that  (^2)  and  (45)  correctly  give  the  order 
of  magnitude  of  Tq.  In  the  theory  of  turbulent  coagulation  of 
aerosols  is  obtained  the  formula  [6] 


(46) 


Possibly,  the  use  for  f(r,  0)  of  a  more  exact  expression  than 
(41)  would  allow  us  to  achieve  better  agreement  of  both  formulas, 

or  even  make  it  possible  to  definitize  (46). 
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When  e  =  10  cm2/s^  and  v  -  0.14  em2/s,  we  find  from  (45) 
t0  =  0.37  s,  which  is  realistic. 

We  showed  that  during  fulfillment  of  condition  (32)  the, 
Chandrasekhar  equation  (25)  for  space-time  correlation  function 
f(r,  t)  has  unique  solution  in  the  class  of  functions  expandable  in 
power  series  in  even  degrees  of  r  and  t,  completely  determined:  by 
assignment  of  space  correlation  function  f(r,  0).  Solution  of 
f  (r,  t)  can  be  found  by  the  formulas  derived  above.  The  form  of 
f(r,  0)  should  be  determined  independently  either  from  theory  or 
from  experiment.  Furthermore,  the  mean  square  of  the  velocity 
component  i?  and  kinematic  viscosity  v  have  to  be  known. 

Thus  the  Chandrasekhar  equation  is  compatible  with  the  well-known 
Kolmogorov  theory,  or  in  general  with  any  theory  giving  f(r,  0),  if, 
of  course,  it  does  not  contain  propositions  contradicting  conditions 
of  homogeneity,  isotrop,  and  stationarity,  or  the  hypothesis  on 
quasi-normality  of  velocity  distribution. 

Certain  theoretical  results  show  that  the  hypothesis  of  quasi¬ 
normality  cannot  be  fulfilled  exactly  (see,  for  example,  work  [4]). 
However,  it  would  hardly  be  correct  to  interpret  such  results  as 
proof  of  the  total  unfoundedness  of  any  derivation  obtained  through 
the  use  of  this  hypothesis,  it  is  known  that  the  system  of  equations 
for  correlation  functions  is  infinite,  so  that  in  any  conceivable 
theory  of  turbulence  it  will  be  necessary  to  adopt  some  or  other 
hypotheses  in  order  to  be  limited  to  a  finite  number  of  equations, 
and  it  is  fully  possible  that  any  of  these  hypotheses  will  lead  to 
more  or  less  considerable  contradictions  and  noncorrespondences. 

4 .  Universal  Solutions  of  the  Chandrasekhar  Equation 

Let  us  examine  a  special  family  of  solutions  of  the  Chandrasekhar 
equation  for  which  we  will  turn  in  equation  (25)  to  dimensionless 
variables 


pSHSaWOTBIfl 


i 


?-+•?**  - 


where  the  coefficient  2  is  introduced  for  considerations  of 
convenience,  mentioned  below,  see  formulas  (102)  and  (103).  Then 
we  have,  assuming  that  f  depends  on  p  and  t. 


0</- 


(48) 


This  equation  does  not  include  quantities  v  and  a*,  characteriz¬ 
ing  liquid  and  turbulence.  Therefore  under  corresponding  boundary 
conditions  it  gives  universal  solutions.  Instead  of  (28)  we  obtain 


-*••**■  -  «,+  ,)<»+.,<*+■>  P  +  »(«  +  *W  +  0  X 
X<4i«  +  32i  +  63K**  +  si0»  +  2)(/»  +  I)(^  +  7)VM«^^],  (49) 

HH 


where  V  —  coefficients  in  expression 


/<*  «>«•  % 


(50) 


If  a  condition  analogous  to  (32)  is  met,  relationships  of  type 
(35)  and  (36)  apply,  where  the  latter  has  the  form 


(51) 


However,  universal  solutions  contradict  the  Kolmogorov  theory 
actually  from  comparison  of  (50)  and  (44)  we  easily  find 


(52) 


where  a„  is  a  dimensionless  quantity,  while  from  the  Kolmogorov 
theory  it  follows  that 


(55) 


where  L  -  magnitude  of  large-scale  vortexes.  Apparently,  universal 
solutions  are  of  little  Interest  from  a  physical  standpoint. 
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5.  Solutions  not  Depending  on  Time  Interval 

Equation  (25)  also  has  solution  not  depending  on  time  interval, 
the  study  of  which  is  o?  certain,  basically  met.  odical,  interest. 

Let  us  record  stationary  equation  (25)  in  the  form 

Representing  f(r)  in  the  form  of  a  series 

M 


(5M 


(55) 


we  find 


fl/+j 


i 

_  2  O'  +  2)(f  + 1)(2^  +  7) 

_ _ _ 

2*«  (i+3){i  +  2)«  +  D(4i*  +  32i+63) 


(56) 


where  i  >_  0.  If  a,  and  a~  are  assigned,  formula  (56)  permits  finding 
all  other  a*  k>Z. 


If  conditions  (32)  is  met,  wltn  the  help  of  formula  (35),  where 
we  consider  «•.*+,  ■•0,  we  car.  determine  a2. 


»*Gi 
280**  ’ 


(57) 


so  that  only  a^  is  arbitrary. 

Equation  ( 44 )  a?^so  has  a  particular  solution  of  special  kind. 

In  order  to  find  it  we  will  start  with  equation  (25),  containing  time. 

Let  us  integrate  both  parts  of  equation  (25)  from  r  to  «  and 
assume  that  (32)  is  satisfied.  We  then  obtain: 

(£-«*)/~  -  i\i?.  (58) 

f 

We  will  represent  solution  of  equation  (58)  in  the  form 
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fir,  0  ~f0(r>  t)  +  u'g(r,  i). 


(59) 


where  fg  satisfies  equation 

•#-- ^o5/,-o. 


with  boundary  condition 

Mr,  0)  -/(*>.  0), 


and  g  satisfies  equation 

- J/(r'.  i)-jpDtn/,t)dr\ 

t 

with  boundary  condition 


g(r,  0)— 0. 

Let  us  find  the  solution  for  (60),  using  apparatus  of  the 
function  theory,  described  in  work  [1],  similarly  to  the  way  in 
this  is  done  in  work  [2].  Let  us  represent  f^  in  the  form 

Mr,t)-/(r,0)  +  f(r,t)0 

where  4  satisfies  equation 

A-*£gf~*££/(r,0) 


with  boundary  condition 


f(r,  0>-0. 


The  solution  of  (55)  is  presented  in  the  form 


r !.  { 


(60) 


(61) 


(62) 


(63) 

S- 

whicn 

(6b) 


(65) 


(66) 


( 67 ) 
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r  yjwg;1 
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Here  G 
space 


—  Green's  function,  dv'  —  element  of  volume  in  five-dimensional 


dxf~  r/*stn*Xsin*1P*ln  Gdr*  dXdVdtdf, 


(68) 


where  r*  -  "distance  to  origin  of  coordinates,"  and  X,  ?,  e  -^nd  <J>  — 
"angles"  formed  fay  radius  vector  with  coordinate  axes  [l]t  Integration 
is  carried  out  over  the  entire  infinite  "space." 


Green's  function  sati.  fies  equation  [1]. 

~-*tAa — a. 


(69) 


using  lyhich,  we  easily  obtain  form  (67) 


t(r.  4 - —Ij/fr-.  0)  Odv'dt'  ~f[r,  0J. 


(70) 


then  from  (64)  we  have 


/•(r,  0-—— rj/(r’,  B)Q(\? —r\, 

j/fr\O)0<t?-rk  flrd&. 


* 


(71) 


The  singular  part  of  Greev-’s  function  is  equal  to 


% 


r^t*,,«**t*«»iieX»ln»Vaine  J 

L  — *-**  j’ 


(72) 


Let  us  assume  that  Green's  function  G  is  different-  from  0  only 
when  t  >  0.  Thez*efore  we  will  assume 


0(r.  9-O,(r.0  +  ~Gt(r.0{2G>^^  J>j>. 


(73) 
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As  a  result  we  find  for  t  >  0 


0(,.o — 


(2*)***/* 


<7*0 


Condition  (73)  together  with  the  condition  of  parity  of  and 
3  *lth  respect  to  t  mean  that  the  sought  solution  of  f(r,  t)  does  not 
have  to  depend  on  t.  Below  we  will  confirm  that  this  is  indeed  the 
case. 


Let  us  represent  f(r,  G)  in  the  form 


/(<',0)-l+a1^+a^+ —  £<V*;«W“I. 


(75) 


The  series  contains  Only  even  degrees  of  r  owing  to  isotropy  [7]  and 
is  assumed  uniformly  convergent  everywhere. 

We  rewrite  (71)  in  the  form 


/•(no— 


i_*_  a 

4»*»! 


,3  5 


0)AT(r',  /)*% 


where  K  is  found  after  integration  with  respect  to  angles 


(76) 


'**-}'**)* 

+-?((''*  *+ J •-"<“) + -  ‘^-T)  - 


(77) 


where 


(76) 


Integrating  in  (76),  we  find 
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0 -j ^  5?’  {«•+  S““[(k  ~ 

~(s+»+'S*)t,K*‘+^'‘l_  >•(»+»>}'  (79> 

where  eQ  —  integration  constant,  proportional  to  p  *,  and  Uk(r,  p)  — 
functions j  equal  to 

U>  «  Jr  *  te-***  +  (- V'''-  ( 80 ) 

Integration  gives 


-2.¥gC^'-^I, 


(81) 


Substituting  (8l)  in  (79),  we  find  an  expression  for  fQ  in  the  form 
of  a  series,  the  structural  rule  for  which  can  oe  seen  if  we  write 
out  several  of  the  first  members: 


Mr.t- i+m+,.)+ai.(-M. 


+  t±j-  +  r 


+  «l 


-( 


£-7-9 


+  3|^-f3 
2-2J* 


2^ 

^+r*U 


’) 


3.7»g»n 

?.2.2.2-^ 


2? 


+  4^|^+6^1  +  4 


2*2*2f* 


2-2^* 


(82) 


The  coefficient  before  the  term  with  ordinal  number  I  in  parentheses, 
ahead  of  which  stands  coefficient  a»>  is  equal  to  the  binomial 
coefficient  C*. 

Let  us  return  to  expression  (59)  for  solution  of  the 
Chendraselchar  equation  (P.5)-  We  will  look  for  g(r,  t)  from  (62)  by 
the  method  of  successive  approximations,  substituting  instead  of 
f(r* ,  t)  a  series  with  unknown  coeffici-ntR 

/(r,  <>-  S  S +  * (?'  r*),  (i3) 

where  »&(<*,  r*)  -  is  an  even  function  of  t  and  r.  Parity  of  function 
?(r,  t)  %»ith  respect  to  t  follows  from  scatlonarity  [?]. 
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Let  us  examine  a  certain  approximation,  for  which 

/(r.  0-  £  £  V*<*:  **-«*• 

MM 

If  coefficients  c#  are  knov/n,  then  g(r,  t)  can  be  recorded  in  the 
form 


(84) 


t  m 


t(r ;  0-jJ[j7(^.  /)  ~~Dt/(f,  /Jrfr'jx 

XO( \?~7\.  |/'-/j)<fe'rfr\ 


(85) 


where  G  —  Green’s  function  (74).  Let  us  replace  variables 


(86) 


Assuming  that  series  (84)  can  be  differentiated  and  integrated 
term-by-term,  we  find: 


J  /{r'  6)  bDtf(r'  *)dr=  C(9)  ~  j  /<r*  °>  ®)dr~ 

r*  0 

-  ■ c«  ■ -  s  s  £  s 


*•  Ji-0  at  .0  «~<1 


(87) 


-*ftere 


^  CO 

C(®)«J/(r.»)~^/(r.  8)</r  = 

/J«0 


(88) 


where  coefficients  C.  are  still  unknown, 


Let  us  find,  auxiliary  function  W(r,  8) 


W(r,  0)-j[j7(r.  0)  JL  DrJ{r%  ft)«/rj;«r',  r)  r'V 


!>  ^ 


O.-J  / 
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Knowing  W(r,  6),  we  can  find  g(r,  t)  by  the  formula 


f 

6{r,t) — ttH:*'*'*  6)rf8‘ 


0  #* 


Substituting  f(r,  t)  here  from  (84),  we  find 


(90) 


Here 


r(r,  »)-i((-2/r)[(6.!«“  +|»Tr*,)c<*>  + 

+*M>»]+j  SSSM- 


k-2  If  0  m~0  K-.V 


P  _2*(a  +  a><»-2) 

*""*  2«+2m-2  X 


wa*+ii»+3 

.i*(2A  +  2«+3) 


+  (a+2*  +  l^r)  £/**«  ~  - 


where  function  U^r,  3)  is  determined  by  (80),  and 

*-v2  1—9  «-0  «-0 

y-  2*(2*  +  3H2fc-2)  ,»+»■  dU*  1 

(2*4-2«-2)(2*  +  2m)  “*** 


(9D 


(92) 


(93) 


It  is  easy  to  verify  that  the  last  expression  equals 

0»—  f {j/<f".  •) p  D./0',  ») **  } r-ir-  -  £  D/*. 

0  >  I*-® 


It  is  possible  to  record  P***  in  the  form 


(94) 


(95) 


where 


7  11  ?,y  Td  <• 


Q-- 


2*(2»4-3H2*~2)  T  2(/n»aw»3 
ai  +  2*-2  [/»<2&  +  2*  + 

1  Z/  1  ' 


■[5fej5i+(^+^t,« 

« — “  4i)  1  ]* 


(96) 


We  will  first  determine  members  of  the  sum  from  (87),  enclosing 
members  with  the  same  highest  degree  of  0  in  parentheses 

£  £  £  £  “  (p»*») + (Pjh» + ^3»i») + 

MMliaO*«4 

+  (^4N0  +  ^JW  +  ^ »M  +  P JIM  +  P JWl)  + 

+  (P$m  +  Pm*  +  P m  +  Pjta*  +  Pnm  +  Pat oi)  +  •••  *“ 

■“  (fhtQn)  +  (<*j» Qt»  +  ®2*^ioQjj)  + 

+  (®«Q«*  +  +  <*»Qa  +  + 

+  (f*mQm  +  awZitQo  +  ^htPaaQn  +  fls*»Q»  +  Qh&Qm  +  ^jo^i^Qj*)  + 

•v 

*»  . 

+  —  ***  k^aaQal  +  ^aao  **"  3joJ  + 

_ * _ 

^  Q<# +  (<*»•+ ^J*0**)  ®7Q»  J  + 

+  (•**!  +  • 

where 


(97) 


+-¥C&*%+'‘> 


P* 


(-A-™.  +  3^-  +3-Z-  +  i'). 

W"  J.  \2.2.2?»  2-2.jP  2?  /’ 


P* 

„  6VT/S.7.9.1I  .  ^7-9‘U*/*  .  *  9-ll  r*  ,  ,  1M  ,  j»\ 

Q* — +4*2*r  +6T2 *T  +*-2T  +  ^) 


(98) 


Let  us  note  that  expressions  in  parentheses  in  (98)  and  (82) 
accurately  coincide,  as  it  should  be  according  to  (79)  ana  (96). 


3? 


We  substitute  (97)  in  (91),  then  substitute  the  obtained 
expression  for  W(r,  9)  in  (90).  After  integration  and  substitution  of 
g(r,  t)  and  fQ(r,  t)  from  (81)  in  (63) 


/  O’,  t)  °)  +  2*/  0  *^10  +  2*7 a*r*  +  3*9*0,/  + 

4*4*lla</  +  5,13,fl«f*  +  ...)  + 

+  +  3*7*9a*r*  +  6*9*  Ua*r4  +  10-11-130*^*  +  ...)  + 

l  \  3*  10  *  *  *  *  ’  / 

+  i*  1 4-  2*70*1-*  4-  3*96/  +  4*116/  +  5*13*6/  + .  * .)  +  0(P) ]. 


From  the  condition  of  stationarity  of  turbulence  it  follows  that 
•odd  degrees  of  t  have  to  drop  from  f(r,  t),  see  (83),  where  this 
requirement  satisfies  each  of  the  approximations  (see  (84)).  Equating 
coefficients  of  members  containing  t  to  zero,  we  find 
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*  3  * 


22-1  + 
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2*3*9* 


2-5.13** 


*  6$< 


4113 


«*- 


a»(2*  +  3)*  1  ^“J,*"^*-"*** 

,  £fie±aiezJ3> 

*’  »(2«  +  3)(«  — 1)  * 


(100) 


It  is  easy  to  verify  that  the  formulas  of  (100),  starting  from 
expression  for  a^Q,  are  equivalent  to  (56),  and  that  the  formula  for 
a2Q  coincides  with  (57). 

If  expressions  for  a^Q  are  substituted  in  (99),  it  will  appear 
that  all  coefficients  on  members  containing  t^,  turn  into  zero,  and 
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thus,  expansion  of  f(r. 


t)  in  series  in  degrees  of  t  has  the  form 


r 

! 

£ 


i 

I 


/(r,  0 -/(r9  0)4-0  (/»). 


(101) 


Examining  (100),  one  readily  sees  that  all  coefficients  akC 
are  proportional  to 


(102) 


Relationships  of  (100)  constitute  recursion  formulas  where  the 
last  coefficient  is  obtained  as  the  sum  preceding  products,  plotted 

in  accordance  with  difined  rule.  It  is  easy  to  see  that  the  factor 

—•  enters  in  expression  for  a^Q  in  degree  k.  Consequently,  series  (75) 
can  ce  rewritten  in  simpler  form  by  introducing  instead  of  r  the 
dimensionless  coordinate 


(103) 


what  corresponds  to  (36).  Then  instead  of  (65)  we  have 


fit*  0)  **  i  4*  «,p*  4-  ' 1  4- 


M 


(104) 


where  coefficients  are  determined  from  recursion  formulas: 


«« 


*■“  .(s(*+3)  ' 

p(2f  +  3Hp—  1) 
«(2«  +  3)(»-*  1)  * 


(105) 
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Expressions  for  6q  and  dg  can  be  found  by  switching  from  r  to 
p  m  formulas  (88)  and  (9*0: 

K —  f{)'m  —DJirt  ( 106 ) 

•  » 


s- 


(107) 


Introducing  dimensionless  time  t  according  to  (47),  we  can 
record  (58)  in  dimensionless  variables: 


(108) 


From  (106),  with  the  help  of  (108),  we  find  a  simpler  expression 
than  (106) 


*. — j|fOS/W«. 


(109) 


where  the  following  relationship  was  used 

|  -o, 


(no) 


following  from  (101). 


With  the  help  of  (110)  it  is  easy  to  show  that  f(p,  t),  in 
general,  does  not  depend  on  t.  Actually,  by  differentiating  (108) 
twice  with  respect  to  t,  we  find 


££\  -m  _ o, 


(in) 


i.e. , 


/<p,*)^/<P.0).  (112) 

Relationships  (105),  (107),  and  (109)  make  it  possible  to  find 
all  coefficients  a,,  inasmuch  as  c  and  n ,  completely  determining  ct,  , 
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are  themselves  determined  by  all  The  value  of  ak  should  decrease 

rapidly  with  increase  of  k,  so  that  it  is  sufficient  to  know  a 
I  certain  quantity  of  first  for  a  given  degree  of  accuracy  in 

I  determination  of  and  f(p,  0). 

j  Coefficients  were  calculated  on  the  "Minsk-l"  computer.  In 

{  integrals  (106.)  and  (109)  instead  of  »  a  certain  was  selected 
!  which  then  was  increased  gradually. 


We  selected  and  arbitrarily  and  then  in  first 

approximation  obtained  etc.,  until  in  a  certain  n-th 

approximation  changes  of  ?  and  n  become  negligibly  small  (with  accuracy 
to  the  fourth  significant  digit).  When  pjn  *  3*5  and  4  process 
converged,  and  when  p  4.5  the  convergence  interval  decreased  so 
much  that  for  5  and  n  found  from  extrapolation  of  curves  and 

(see  Fig.  1)  divergence  of  the  iterative  process  was  observed  so  that 
it  was  necessary  to  select  such  a  and  pair  for  which 

divergence  was  the  slowest.  Figure  1  shows  how  stabilization  of 
C  and  n  occurred  with  growth  of  pm. 


The  most  probable  value  of  C  and  n  are  Q  =  0.3585  and 
n  =  -0.04154,  with  accuracy  to  2-3  units  of  the  fourth  digit.  With 
growth  of  p  the  number  of  N  of  coefficients  necessary 
for  finding  z,  and  r\  with  an  accuracy  to  one  unit  of  the' 
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fourth  digit  increases  sharply, 


With  =  5.95  are  N  =  700  was 
m 


required.  Further  increase  of  pm  turned  out  ..to  be  impossible  in 
connection  with  the  limited  capacity  of  the  fast  store  of  the 
machine  (2048  numbers)  and  increase  of  computation  time.  For 
obtaining  ?  and  n  from  values  of  the  preceding  approximation  two 
hours  of  continuous  machine  operation  were  required. 

From  the  found  final  values  of  £  =  0.3585  and  n  «  -0.04154 

a  graph  of  f(p,  0)  was  plotted,  see  Fig.  2.  From  the  graph  It  is 

clear  that  f(p,  0)  intersects  the  axis  of  abscissas  when  p  'v  6. 

Values  of  f(p,  0)  for  c  =  0.3585  and  n  *  -0.04157  differ  from  the 

above-mentioned  by  less  than  on  0.001  with  p  _<  6,  which  permits 

judging  the  accuracy  of  plotting  of  the  graph.  In  Fig.  2  is 

2 

also  given  a  graph  of  function  expC-ct^p  )  for  =  0.07170.  It  is 
clear  that  up  to  p  =  5  both  curves  almost  coincide,  so  that  decrease 
of  f(p,  0)  turns  out  to  be  very  rapid.  The  found  solution  is  checked 
by  direct  substitution  of  f(p,  0)  in  the  form  of  a  series  in  equation 
(108),  for  p  *  1. 


Fig.  2. 


9  2*3* 


The  obtained  solution  for  f(p,  0)  is  interesting  basically  from 
methodical  stand  point.  Most  interesting  is  the  qult-e  rapid 
convergence  of  the  series  representing  the  solution,  which  permits 
relying  on  analogous  convergence  of  series  giving  both  time-independent 
and  time-dependent  solutions.  Green's  function  for  the  linear  part 


of  the  Chandrasekhar  equation  (25),  used  during  the  finding  of  f(p,  0) 
can  be  useful  during  more  detailed  investigation  of  this  equation. 

6.  Conclusions 

The  conducted  investigation  lead  to  the  following  results. 

:.)  A  dynamic  equation  was  derived  for  determining  of  the 
scalar  of  correlation  tensors  describing  two-point  space-time 
correlation  of  velocities  of  turbulent  flow  and  being  a  generalization 
of  the  Chandrasekhar  equation  for  stationary  turbulence.  During 
derivation  we  adopted  conditions  of  homogeneity  and  isotropy,  and 
also  the  generalized  hypothesis  on  quasi-normality  of  velocity 
distribution. 

2)  During  direct  derivation  of  the  stationary  equation  a  pair 
of  equations  not  noticed  by  Chandrasekhar  is  obtained.  However,  these 
equations  one  should  be  rejected,  inasmuch  as  it  is  impossible  to 
obtain  them  from  the  more  general  nonstationary  equation. 

3)  In  principle  it  is  possible  to  obtain  one  more  equation  for 
nonstationary  turbulence,  which  has  not  yet  been  investigated. 


*0  The  generalized  Chandrasekhar  equation  correctly  describes 
the  last  stage  of  degeneration  of  turbulence. 

The  Chonri-asekhar  equation  (stationary)  for  the  space-time 
U  if  .git  . . -onl  correlation  function  can  be  solved  by  way  of  the 

..  .-ntation  of  solution  in  the  form  of  a  power  series.  The  solution 
•Is  unique  if  one  assigns  purely  spatial  ana  tcmporarl  correlation 
functions,  which  can  be  determined  by  another*  theory  by  experiment. 


6)  If  a  certain  known  combination  of  partial  derivatives  of 
t:.«  ..orrelation  function  turns  to  zero  at  infinity,  for  unique 
solution  of  the  Chandrasekhar  equation  it  is  sufficient  to  assign  tr. 
j_...  correlation  function.  Here  the  Euler  '  tae  tempera :  rlc- o.::-  :  . 


of  turbulence  turns  out  to  be  a  function  of  the  mean  square  of 
velocity,  kinematic  vi  cosity,  and  the  integral  of  known  form  from  the 
space  correlation  function. 

7)  The  Chandrasekhar  equation  has  a  class  of  solutions,  the 
form  of  vjhich  is  universal,  i.e.,  does  not  depend  on  the  properties  of 
liquid  and  intensity  of  turbulence.  These  solutions  do  not  agree 
with  the  formula  for  the  mean  square  of  velocity  of  turbulent  flow, 
known  from  the  Kolmogorov  theory. 

8)  We  investigated  solutions  of  the  stationary  equation 
not  depending  on  time  interval,  and  in  particular,  one  singular 
solution,  determined  with  the  help  of  Green’s  function. 

In  conclusion  we  must  point  out  that  results  of  ’‘his  work 
can  be  used  during  derivation  and  solution  of  equations  for  other 
space-time  correlation  functions  in  turbulent  flow. 

The  authors  thank  A.  I.  Ivanovskiy  for  useful  criticism  and 
Yu.  V.  Sidorov  for  discussion  of  a  number  of  questions  of  mathematical 
character. 
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TURBULENCE  SPECTRUM  OP  A  STABLY 
STRATIFIED  ATMOSPHERE 


\ 


\ 

\ 


G.  .  .  Shut* 


Results  of  experimental  investigations  ol 
the  turbulence  spectrum  in  the  stratosphere  and 
upper  troposphere#  where  temperature  stratfi- 
catlon  is  stable,  are  examined  a  diagram  of 
the  generalized  turbulence  spectrum,  including 
the  buoyancy  subdomain  is  offered.  Particular 
cases  of  spectra  with  limited  band  of  wave  numbers, 
including  the  spectrum  of  nquasi-waven  distur¬ 
bances  are  examined.  For  a  physical  interpre¬ 
tation  of  energy  transitions  in  the  spectrum 
curves  «9~/(a)  are  used. 


I.  Introduction 


Turbulence  in  the  free  atmosphere  is  the  phenomenon  determining 
in  many  respects  the  physics  of  processes  occuring  in  it.  Furthermore, 
study  of  atmospheric  turbulence  is  necessary  for  solution  of  a  whole 
series  of  purely  applied  problems  of  dynamics  of  flight  and  aircraft, 
of  propagation  and  scattering  of  purities  in  the  atmosphere,  etc. 

Free  atmosphere  is  usually  stratified  stably  in  density,  which 
noticeably  affects  the  character  of  appearance  and  development  of 
turbulence  in  it.  The  question  of  appearance  of  turbulence  in  a 
thermally  stratified  medium  was  solved  theoretically  by  Richardson, 
who  gave  the  criterion  of  growth  and  decrease  of  turbulent  energy : 
turbulence  increases  when  R^  <  1  it  is  not  changed  when  R.  =1,  and  it 
fades  when  R.  >  1,  where 
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(1) 
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T«~T 


Here  from  equation  (X)  it  is  clear  that  Richardson  considered 
both  thermal  and  dynamic  stratification.  However,  for  steady 
turbulence,  i.e.,  the  Richardson  number,  does  not  permit  estimating 
the  Intensity  of  turbulence  and,  this  is  especially  important,  gives 
no  information  about  scales  of  turbulent  motion. 


Since  turbulence  is  characterised  by  field  of  random  veloeiting, 
for  description  of  the  field  of  turbulence  statistical  characteristics 
are  used.  One  of  such  character!?  ies,  allowing  us  to  obtain 
information  on  the  intensity  of  t.  ">ulent  motion  of  different  scales, 
is  the  spectral  density  of  energy  distribution  of  turbulent  motion, 
or  the  energy  spectrum  of  turbulence 


5(C)  *- 


\ 

2  *8* 


(2) 


where  jj  I  -  scale  of  motion,  -  vorticity  of  corresponding  scale. 

For  locally- isotropic  uniform  turbulence  Kolmogorov  and 
Obukhov,  proceeding  from  considerations  of  self-similarity  and  assuming 
that  energy  is  transmitted  from  larger  scales  to  smaller  without 
losses,  obtained  an  expression,  which  is  spectral  form  is  known  as 
the  ”-5/ 3"  law 

.1  ± 

5(C) w«3S  \  (3) 


where  e  is  the  rate  of  energy  transfer  over  the  spectra,  numberically 
equal  to  the  rate  of  dissipation  of  turbulent  energy  as  heat.  As  car. 
be  seen  from  expression  (3),  e  also  characterizes  the  intensity  of 
turbulence . 


We  will  assume  that  ir:  turbulent  atmospnere,  even  if  we  place 
no  limitations  on  spatial  structure  and  distribution  of  turbulent 
formation  (i.e.,  also  in  the  case  of  anisotropic  and  nonuniform 
turbulence),  energy  transfer  from  larger  scales  to  smaller  takes 
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place.  In  the  language  of  spectral  concepts  this  means  that  energy 

flux  in  the  spectrum  is  directed  toward  greater  wave  numbers.  Hence 

follows  directly  determination  of  the  rate  of  energy  transfer  in  a 

spectrum  i,as  a  function  of  wave  number  q  .  In  particlar  case  (3) 

const.  Here,  as  already  indicated,  in  a  certain  interval  of 

wave  numbers  energy  does  not  enter  and  is  not  expended,  but  is  only 

transferred  from  some  movements  to  others.  In  Pig.  la  in  coordinates 

S(Q)  and  fi  are  depicted  e  =  const  lines,  where  e,  <  s~  <  e0.  If  the 

JL  c  j 

spectrum  of  turbulence  is  in  parallel  with  these  lines,  it  means  that 
the '”-5/3”  law  (segment  BC)  is  valid.  If,  however,  the  spectral 
curve  with  increase  of  wave  number  passes  from  one  e  to  another 
(segment  AB),  it  means  that  in  a  certain  interval  of  wave  numbers 
energy  enters  spectrum  and,  conversely,  in  the  interval  corresponding 
to  segment  CD  kinetic  energy  of  turbulent  pulsations  does  not  pass 
completely  into  energy  of  pulsations  of  smaller  scales.  In  other 
words,  in  region  A^B1  the  source  acts,  but  in  region  the  consumer 

of  energy  of  turbulent  pulsations  is  active. 


Fig.  1.  Spectral  density  of  the 
energy  of  turbulent  pulsations  and 
the  rate  of  energy  transfer  over 
the  spectrum. 

A  very  convenient  form  of  representation  of  experimental  data  in 
those  cases  when  we  want  to  trace  transfer  of  energy  in  the  turbulence 
spectrum  is  a  graph  of  function t#=/ (G),  presented  in  Pig.  lb.  On 
such  a  graph  it  is  possible  to  separate  distinctly  both  regions  of 
entry  of  energy  into  spectrum  and  regions  where  enegry  is  "sucked” 
from  the  spectrum. 
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II .  Turbulence  Developed  in  Stably 
Stratified  Atmosphere 


Stable  stratification  of  medium  prevents  the  appearance  in  it 
of  oscillations  whatever.  If,  however,  in  such  medium  chaotic 
turbulent  movements  nevertheless  appear,  tnen,  according  to  the 
Richardson  number,  this  means  that  the  dynamic  factor  (in  the  opinion 
of  Richardson  this  factor  is  the  vertical  velocity  gradient) 
predominates  over  thermal  stability.  However,  the  influence  of 
Archimedian  forces  is  not  limited  by  the  fact  that  they  prevent  the 
the  appearance  of  turbulence.  In  well-developed  turbulent  flow 
Archimedian  forces  influence  the  character  of  the  turbulence  spectrum. 
The  range  of  wave  numbers  (scales)  in  which  this  influence  is 
significant  is  called  the  subdomain  of  buoyancy.  Vortices  in  tn  s 
region  during  their  lifetime  must  accomplish  work  against  the 
Archimedian  forces,  on  which  they  expend  part  of  their  kinetic  energy, 
i.e.,  forces  of  negative  buoyancy  are  the  consumer  of  turbulent 
energy  here.  As  was  indicated  above,  *a  drops  with  increase  of  2, 
and  the  spectral  curve  has  slope  greater  than  5/3. 

First  experimental  data  confirming  the  presence  in  the 
turbulence  spectrum  of  an  interval  of  wave  numbers  in  which  the  forces 
6?  ■buoyancy  appear  were  obtained  in  1959-1960  [3],  [4], 

Analyzing  results  of  radio  experiments,  Bolgiano  [4]  offered 
a  certain  theoretical  model  of  the  spectrum  of  turbulence  ir.  the 
buoyancy  interval.  Assuming  that  in  a  defined  range  of  wave  numbers 
significant  influence  on  the  form  of  the  spectrum  is  rendered  by  the 
rate  of  dissipation  of  mean  square  fluctuations  of  specific  forces  of 
buoyancy,  Bolgiano  drew  a  conclusion  to  the  effect  that  in  the 
interval  of  buoyancy  the  slope  of  the  spectral  curve  remains  constant 
but  different  from  5/3-  He  obtained  expressions  for  spectra  of 
pulsations  of  velocity  and  temperature 

_n 

Sr(«)cc2~'‘.  (4) 

T_ 

ST<S)ccSf‘,  (5) 
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i  where  S.(Q)  —  energy  spectrum  of  velocities  of  turbulent  pulsations, 

i  $,(&)  -  energy  spectrum  of  pulsations  of  temperature.  During  the 

i 

analysis  of  results  of  experimental  investigations  of  the  spectrum  of 
the  vertical  component  of  turbulent  pulsations  of  velocity,  which  were 
conducted  at  TsAO  in  1959-1960,  it  was  found  that  in  the  interval  of 
wave  numbers  corresponding  to  scales  of  from  hundreds  of  meters  to  two 
or  three  kilometers  the  experimental  curve  has  slope  considerably 
greater  than  5/3  (see  article  [3]). 

Investigations  were  made  with  help  of  a  flying  laboratory  in 
the  upper  troposphere  basically  in  clear  sky.  Zones  of  intense 
turbulence  in  jet  streams  were  inspected.  Thermal  stratification  of 
the  atmosphere  at  these  heights  was  stable. 

In  work  [3]  an  expression  was  offered  for  the  turbulence 
spectrum,  agreeing  well  with  experimental  data 

JL  A  _± 

S.(8)co*#32~  a(i  +  62  »),  (6) 

where  Sa{Q)  —  energy  spectrum  of  vertical  pulsations  velocity,  eQ  ~ 
rate  of  dissipation  of  turbulent  energy  as  heat,  b  -  coefficient, 
depending  on  gradient  of  potential  temperature,  i.e.,  on  the  degree 
of  thermal  stability  of  the  atmosphere. 

r 

i 

Expression  (6)  was  obtained  from  the  assumption  that  the  rate  of 
conversion  of  energy  in  the  subdomain  of  buoyancy  depends  on  the  wave 
<  .  number  and  is  determined  by  the  average  gradient  of  buoyancy  forces. 

As  was  shown  in  the  cited  work,  expression  (6)  is  not  strict.  The 
problem  of  the  turbulence  spectrum  was  more  strictly  theoretically 
solved  by  Lumley  [6].  He  started  from  the  same  considerations  on 
dependece  of  the  rate  of  transfer  of  energy  in  the  spectrum  on  the 
gradient  of  buoyancy  forces,  as  in  article  [3],  and  obtained  the 
expression: 

J2_  _  * 

L 
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(7) 


where  Q*  is  wave  number,  characterizing  the  subdomain  of  buoyancy. 
Expression  (7)  coincides  completely  with  expression  (6).  Thus  it  is 
possible  to  consider  experimentally  and  theoretically  established  the 
presence  in  the  spectrum  of  developed  turbulence  during  stable 
stratification  of  the  atmosphere  of  a  subdomain  of  buoyancy,  which  is 
characterized  by.  the  slope  of  the  spectral  curve  greater  than  5/3* 

The  distinction  between  expression  (4)  on  the  one  hand  and  (6) 
and  (7)  on  the  other  is  reduced,  as  a  rule,  to  distinction  in  the  value 
of  the  exponent  on  ft  in  the  subdomain  of  buoyancy ,  so  that  in  the  first 
case  the  slope  of  the  spectral  curve  is  constant,  while  In  the  second 
it  is  a  function  of  wave  number. 


In  recent  years  in  TsAO  much  experimental  data,  obtained  under 
conditions  when  it  was  possible  to  expect  the  presence  of  the  subdomain 
of  buoyancy  in  the  turbulence  spectrum,  has  been  accumulated.  In 
Fig.  2  are  presented  curve  of  e^,  obtained  in  1565  during  flights 
the  flying  laboratory  in  clear  sky.  Curves  1  and  2  were  obtained 
for  the  horizontal  component.  Horizontal  fluctuations  flow  velocity 
were  measured  in  the  region  of  large  scales  by  a  Doppler  system  [2] 
and  in  the  region  of  small  scales  by  an  aircraft  hot-wire  anemometer 

Cl]. 


Fig.  2.  Experimental  curves 
«a ~/<c>;  1  and  2  —  for  horizontal 
and  3  —  for  vertical  components. 


Curve  3  was  obtained  for  the  vertical  component  according  to 
measurements  of  vertical  overload  of  the  center  of  gravity  of  the 
aircraft,  with  account  being  taken  of  its  transfer  function.  The 
curve  liec  considerably  higher  than  1  and  2,  since  it  corresponds  to 
the  presence  of  intense  bumping  of  the  aircraft,  i.e.,  to  considerably 
higher  energy  of  turbulent  motion. 


The  curves  in  Pig.  2  were  obtained  during  stable  temperature 
stratification,  and  on  them  is  seen  distinctly  a  region  of  decrease 
of  (subdomain  of  buoyancy).  It  is  interesting  to  note  that  for 
greater  intensities  of  turbulence  this  region  shifts  in  the  direction 
of  greater  wave  numbers. 


Experimental  investigations  of  atmospheric  turbulence  widly 
developed  during  recent  years  made  it  possible  to  accumulate  a  large 
quantity  of  data  on  the  structure  of  the  field  of  turbulence.  It  is 
obvious  that  every  method  permits  obtaining  the  spectrum  of  turbulence 
in  only  a  defined  interval  of  scales,  but  at  present  our  concep  ,  of 
the  spectrum  of  turbulence  of  free  atmosphere  frost  scales  equal  to 
tens  of  kilometers  to  submolecular  scales,  where  kinetic  energy  of 
turbulent  pulsations  passes  to  kinetic  energy  of  molecules.,  i.e-.,  to 
heat,  has  become  complicated.  In  Fig.  3  is  represented,  schematically 
the  energy  spectrum  of  turbulence  of  stably  stratified  atmosphere. 


Pig.  3*  Generalized, 
curves  of  s<a>  and 

j.C 

for  stably  stratified 
atmosphere . 


Region  I  is  the  region  of  scales  in  which  primary  turbulent 
formations  appear.  Average  motion  in  these  scales  loses  its  stability, 
and  part  of  the  energy  of  basic  flow  goes  to  formation  of  disordered 
fluctuations.  Region  I  is  characterized  by  growth  of  efi.  This  does 
not  mean  that  only  this  region  receives  turbulent  energy.  The 
increase  of  makes  it  possible  to  conclude  only  that  entry  of 
energy  into  the  spectrum  is  faster  than  the  draining  off  (consumption) 
of  this  energy.  We  will  not  stop  here  on  consideration  of  possible 
mechanisms  of  generation  of  turbulent  energy  in  region  I.  This 
question  is  very  complex  and  is  still  far  from  being  fully  studied. 

Region  II  is  characterized  by  quasi-equilibrium  between  entry 
of  energy  into  the  spectrum  and  loss  of  this  energy.  As  can  be  seen 
from  the  figure,  the  rate  of  energy  transfer  over  the  spectum  e~ 
in  this  region  remains  almost  constant. 

The  following  region  III  —  subdomain  of  buoyancy  —  is 
characterized  by  the  fact  that  at  wave  numbers  2,-2,  acts  a  powerful 
consumer  of  turbulent  energy  while  entry  of  energy  into  the  spectrum 
from  without  on  these  scales  is  practically  lacking.  Owing  to 
loss  of  energy  on  work  against  forces  of  buoyancy,  the  rate  of  energy 
transfer  over  the  spectrum  in  this  subdomain  drops. 

Region  IV  Is  the  classical  inertial  interval.  Here  there  is 
neither  influx  or  loss  of  energy.  The  rate  of  transfer  of  energy 
over  the  spectrum  in  inertial  interval  is  constant  and  is  equal  to 
the  rate  at  which  turbulent  energy  passes  over  the  right  boundary  of 
the  inertial  to  thermal  intervals  —  rate  of  dissipation  eQ. 

Region  V  is  the  viscous  interval,  where  kinetic  energy,  the 
energy  of  turbulence,  is  converted  to  heat.  In  this  region  the  rate 
of  energy  transfer  over  the  spectrum  drops  to  zero. 

In  accordance  with  changes  of  rate  of  transfer  of  energy  over 
the  spectrum  in  different  regions  the  form  of  spectral  curve  will  also 
be  changed. 
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From  the  assumption  about  homogeneity  and  local  isotropy 
Kolmagorov  and  Obukhov  obtained  an  expression  for  the  energy  spectrum 
of  turbulence  in  the  inertial  interval  -  the  "-5/3"  spectral  law 

2.  JL 

S(8)«o«03e~  *.  (8) 

If  by  definition  in  the  inertial  interval  there  are  no 
sources  and  consumers  of  turblent  energy,  but  only  inertial  transfer 
of  energy  over  the  spectrum  from  larger  scales  to  smaller  takes  place, 
then  in  the  interval  of  equilibrium  (region  II)  both  sources  and 
consumers  of  turbulent  energy,  which,  however,  compensate  one  another, 
are  at  work.  In  spite  of  the  fact  that  in  the  interval  of  quasi¬ 
equilibrium  the  condition  of  local  isotropy  is  known  not  to  be  met, 
the  spectrum  in  this  interval  also  is  described  by  the  expression: 

S(8)coC~  *.  (9) 

The  expression  for  the  spectrum  in  the  viscous  interval  is 
obtained  from  the  Geisenberg  equation  [5]  and  has  the  form: 

smooQr*.  (10) 

As  regards  the  subdomain  of  buoyancy,  the  spectral  curve  here  has  slope 
greater  than  5/3  and  is  described,  according  to  [3],  by  expression: 

*  1 

S(0)co2~T(i  +  tff J).  (U) 

III.  Discrete  Spectra.  Wave  Disturbances  in 

the  Atmosphere 

Everything  said  above  pertains  to  continuous  spectra.  However, 
among  experimentally  obtained  spectra  there  also  are  such  in  which 
turbulent  motion  in  a  wide  range  of  scales  is  lacking.  More  precisely, 
this  motion  have  amplitudes  smaller  chan  the  threshold  of  sensitivity 
of  the  measuring  equipment.  Such  discrete  spectra  are  due  to  the 
presence  of  stable  wave  disturbances,  not  transmitting  their  energy 
to  smaller  scales.  s 
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Sometimes  the  spectrum,  continuous  at  wave  numbers  2<2*,  3harply 
drops  when  C>0#.  In  this  case  one  should  either  recognize  that  the 
assumption  about  energy  transfer  in  the  continous  spectrum  from  larger 
scales  to  smaller  is  not  satisfied  or  should  assume  the,  presence  of 
conversion  of  kinetic  energy  pulsations  when  &>fV  t0  som®  other 
form  of  energy. 


We  will,  as  before,  consider  that  in  turbulent  flow  energy 
transfer  from  larger  scales  to  smaller  takes  place.  As  already 
stated,  growth  of  the  rate  of  energy  transfer  indicates  the 
presence  of  a  source  of  turbulent  energy  and,  conversely,  a  decrease 
of  indicates  the  presence  of  a  consumer. 

When  spectrum  of  turbulence,  and  consequently  and  *,>•■» f(S),  are 
continuous,  such  description,  in  general,  is  trivial.  However,  quite 
another  picture  is  obtained  if  with  the  same  assumption  about  the 
direction  of  energy  flow  in  the  spectrum  we  turn  to  consideration  of 
stable  wave  disturbances  in  a  real  medium,  i.e.,  in  a  medium 
possessing  dissipative  properties  (for  example,  viscosity). 


I 


Let  us  examine  forced  oscillations  of  air  with  frequency 
corresponding  to  wave  number  ft.  We  know  that  if  the  source  of  these 
oscillations  ceases  to  act,  the  oscillations  will  fade  with  time 
owing  to  internal  friction  (viscosity)  of  air.  Let  us  try  to  analyse 
the:  meaning  of  this. 

An  air  mass  has  defined  viscosity,  depending  on  temperature  and 
density.  Molecules  of  air  are  in  random  motion,  where  the  number  of 
molecules  in  a  unit  -volume  and  their  average  velocity  determine  the 
viscosity  of  the  air. 


;I  When  we  say  that  oscillation  of  air  fade,  it  means  that  the 

mechanical  energy  of  these  oscillations  turns  to  energy  of  molecular 
motion. 

If  continuous  wave  oscillations  exist  in  the  atmosphere,  it 
means  that  they  are  supported  constantly  by  an  active  external  source. 

5Q 
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The  scale  of  disturbances  generated  by  this  source  should ,  in  any  case;, 
be  no  less  than  the  length  of  the  continuous  wave.  Consequently,  in 
case  of  stable  wave  disturbances  also  the  spectra!  energy  flow 
preserves  its  direction  toward  larger  wave  numbers,  and  energy  going 
to  larger  scales  is  transferred  to  the  region  of  submolecular 
scales.  The  spectrum  in  this  case  is  discrete  and  has  the  form 
presented  in  Pig,  k. 


Pig.  M.  Spectrum  of  purely 
wave  disturbance  in  viscous 
medium. 


On  the  axis  of  abscissas  to  the  first  peak  corresponds  the 
wavelength  of  stable  oscillation;  to  the  second  corresponds  the 
value  of  the  order  of  the  mean  free  path  of  molecules. 

Thus  if  the  medium  possesses  viscosity,  "sucking  out"  of 
energy  of  any  mechanical  motion  appearing  in  this  medium  will  occur. 
It  is  possible  to  apply  the  same  reasoning  to  a  medium  in  which 
random  motion  with  scales  much  larger  than  molecular  exists,  that  is, 
to  a  medium  possessing  turbulent  viscosity.  If  in  an  air  mass  there 
already  are  movements  of  defined  scales,  the  intensity  of  these  move¬ 
ments  can  increase  directly  as  a  result  of  energy  cf  other  mechanical 
motion,  even  if  its  scale  is  much  greater.  Prom  this  stand  point 
it  is  easy  to  explain  the  character  of  the  spectral  curve  obtained 
by  Van  der  Hoven  [73,  having  dips  over  a  wide  range  of  wave  numbers  * 

Let  us  now  consider  cases  when  the  3peetrum  is  continuous  over 
a  rather  wide  range  of  wave  numbers,  but  then  drops  sharply. 

In  examining  of  such  rapidly  dropping  spectra  we  will,  as 
earlier,  consider  that  spectral  energy  flow  is  directed  toward  larger 
wave  numbers..  Inasmuch  as  such  spectra  are  characteristic  for  stably 
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stratified  atmosphere,  one  should  assume  that  rapid  drop  of  spectrum 
is  explained  by  the  influence  of  forces  of  buoyancy,  which  consume 
almost  all  the  energy  arriving  from  the  direction  of  smaller  wave 
numbers. 


IV.  Turbulence  in  a  Thermally  Stable 
'Stratified  MetSCum 


With  strong  thermal  stability  of  atmosphere  which  for  example, 
is  typical  of  the  stratosphere,  for  disturbance  of  Xaroinarity  flow 
large  wind  gradients  are  necessary.  The  ivary  fact  of  appearance  and 
development,  of  turbulence  in  such  a,  medium  indicates  that  this 
condition  Cl)  is  fulfilled. 

However,  although  the  reserve  of  thermal  stability  it 

insufficiently  to  impede  generation  of  turbulence,  the  forces  <->f 

».  ’*  .  * 

buoyancy  intensively  counteract  cascade  energy  transfer  over  the 
spectrum.  Three  cases,  corresponding  to  the  three  curves  in  Pig.  5 
are  possible  here. 


Pig,  5*  Modifications  of 
turbulence  spectra  for 
different  relationships 
between  intensities  of 
source  and  consumer  of 
turbulent  energy  and  for 
their  different  character¬ 
istic  scales. 


li  With  very  strong,  thermal  stability  of  medium  the  spectrum 
of  turbulence  can  be  localized  in  a  narrow  band  of  wave  numbers.  In 
this  case  the  spectrum  obtained  as  a  result  of  treatment  of  realization 
of  finite  length  for  practical  purposes  will  not  differ  from  the 
spectrum  of  purely  harmonic  oscillation,  if  we  use  existing  methods 
of  statistical  computer  treatment.  Such  '“quasi-wave"  disturbances 
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differ  from  those  of  "pure  wave”  type  by  the  fact  that  in  them  energy 
transfer  from  larger  scales  to  smaller  takes  place,  and  then  all 
energy  is  expended  on  work  against  the  forces  of  buoyancy. 

2.  With  somewhat  lesser  thermal  stability  of  atmosphere,  as 
well  as  with  its  greater  dynamic  instability,  and  also  in  those  cases 
when  the  source  of  turbulent  energy  is  located  in  the  region  of  wave 
numbers  distant  from  the-  subdomain  of  buoyancy,  the  spectrum  of 
turbulence  is  found  considerably  broader  and  differs  significantly 
from  the  spectrum  of  purely  harmonic  oscillations.  However,  even  in 
this  case  the  forces  of  buoyancy  absorb  all  the  kinetic  energy  of 
turbulent  pulsations. 

3.  The  most  frequently  encountered  case  is  that  when  there  is 
a  continuous  spectrum  of  turbulence,  in  which  kinetic  energy  is 
transferred  from  larger  scale  to  smaller.  Part  of  the  energy  becomes 
potentia?.  energy,  and  then  in  scales  where  the  forces  of  buoyancy  no 
longer  play  an  essential  role  there  is  inertial  transfer  in  accordance 
with  the  "“5/3*’  law.  In  the  end  all  energy  turns  to  heat  in  the 
viscous  interval. 

Diverse  variants  of  spectra,  shown  in  Fig.  5,  are  particular 
cases  of  the  general  diagram,  presented  in  Fig.  3,  and  characterize, 
distribution  of  energy  in  the  spectrum  of  developed  turbulence,  i.e., 
the  stationary  case.  During  .experimental  investigations  of  turbulence 
it  sometimes  is  possible  to  obtain  spectra  of  developing  or  fading 
turbulence.  Such  spectra  have  complicated  form,  and  on  them  are 
distinctly  seen  local  (according  to  wave  numbers)'  sources  of 
turbulent  energy,  such,  for  example,  as  dis intigrating  gravitational 
waves. 


For  the  spectrum  ,of  developed  turublence  in  a  stably  stratified 
medium  it  is  characteristic  that  the  source  of  turbulent  energy 
lies  in  the  region  of  small  wave  numbers,  while  the  spectral  curve 
itself  has  only  one  peak  in  the  region  of  generation. 
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